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Abstract 

The Tutte polynomial of a graph G is a two-variable polynomial T(C7; x, y) that 
encodes many interesting properties of the graph. We study the complexity of the 
following problem, for rationals x and y: given as input a planar graph G, determine 
T(G; x,y). Vertigan completely mapped the complexity of exactly computing the 
Tutte polynomial of a planar graph. He showed that the problem can be solved 
in polynomial time if (x,y) is on the hyperbola H q given by (x — l)(y — 1) = q 
for q = 1 or q = 2 or if (x,y) is one of the two special points (x,y) = (—1,-1) 
or (x,y) = (1,1). Otherwise, the problem is #P-hard. In this paper, we consider 
the problem of approximating T(G;x,y), in the usual sense of "fully polynomial 
randomised approximation scheme" or FPRAS. Roughly speaking, an FPRAS is 
required to produce, in polynomial time and with high probability, an answer that 
has small relative error. Assuming that NP is different from RP, we show that there 
is no FPRAS for the Tutte polynomial in a large portion of the (x, y) plane. In 
particular, there is no FPRAS if x > 1, y < — 1 or if y > 1, x < — 1 or if x < 0, y < 
and q > 5. Also, there is no FPRAS if x < 1, y < 1 and q = 3. For q > 5, our result 
is intriguing because it shows that there is no FPRAS at (x, y) = (1 — qj (1 + e), —e) 
for any positive e but it leaves open the limit point e = 0, which corresponds to 
approximately counting (/-colourings of a planar graph. 

*This paper is available on the ArXiv at |http: / /arxiv.o rg/ abs / 907. 1724| The work was partially 
supported by the EPSRC grant The Complexity of Counting in Constraint Satisfaction Problems 
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1 Introduction 



1.1 The Tutte Polynomial 

The Tutte polynomial of a graph G = (V, E) (see [TTl 13J) is the two- variable polynomial 

T(G; Xj y)=J2(x- lf V ^~^ E \y - (1) 
ACE 

where k(V, A) denotes the number of connected components of the graph (V, A) and n = 
\V\. Following the usual convention for the Tutte polynomial [TD] a graph is allowed to 
have loops and/or multiple edges. 

Many interesting properties of a graph correspond to evaluations of the Tutte polyno- 
mial at different points (x,y). For example, the number of spanning trees of a connected 
graph G is T(G; 1,1), the number of acyclic orientations is T(G; 2, 0), and the reliability 
probability R(G; p) of the graph is an easily-computed multiple of T(G; 1,1/(1— p)). For a 
positive integer q, the Tutte polynomial along the hyperbola H q given by (x — 1) (y — 1) = q 
corresponds to the Partition function of the g-state Potts model. See Welsh's book [13] for 
details. 

Two particularly interesting Tutte invariants correspond to evaluations along the x axis 
and the y axis. In particular, 

• The chromatic polynomial P(G; A) of a graph G with n vertices, m edges and k 
connected components is given by 

P(G; A) = (-l) n - k X k T(G; 1 - A, 0). 

When A is a positive integer, P(G; A) counts the proper A-colourings of G. 

• The flow polynomial F(G; A) is given by 

F(G; A) = (-l) m - n+fc T(G; 0, 1 - A). 
When A is a positive integer, F(G; A) counts the nowhere-zero A-flows of G. 

1.2 Evaluating the Tutte Polynomial 

For fixed rational numbers x and y, consider the following computational problem. 

Name. Tutte(x, y). 
Instance. A graph G = (V, E). 
Output. T(G] x, y). 

The parameters x and y are fixed in advance and are not considered part of the problem 
instance. Each choice for x and y defines a distinct computational problem. 
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Jaeger, Vertigan and Welsh [6j have completely mapped the complexity of Tutte(x, y). 
They have shown that Tutte(x, y) is in FP for any point (x, y) on the hyperbola Hi and 
when (x, y) is one of the special points (1, 1), (0, —1), (—1, 0), and (—1, —1). They showed 
that Tutte(x,?/) is #P-hard for every other pair of rationals {x,y). See [H] for definitions 
of FP and jfP; informally, FP is the extension of the class P from predicates to more 
general functions, and #P is the counting analogue of NP. Jaeger et al. also investigated 
the complexity of evaluating the Tutte polynomial when x and y are real or complex 
numbers, but that is beyond the scope of this paper. 

Vertigan [12] considered the restriction of Tutte(o;, y) in which the input is restricted 
to be a planar graph. 

Name. PlanarTutte(x, y). 
Instance. A planar graph G = (V, E). 
Output. T{G; x, y). 

He showed that PlanarTutte(x, y) is in FP for any point (x,y) on the hyperbolas Hi 
or H 2 , and when (x,y) is one of the special points (1,1) and (—1,-1). He showed that 
PlanarTutte(x, y) is #P-hard for every other pair of rationals (x, y). The hyperbola H 2 
is of particular interest, as the Tutte polynomial here corresponds to the partition function 
of the celebrated Ising model in statistical physics. 

1.3 Approximating the Tutte polynomial 

A fully polynomial randomised approximation scheme (FPRAS) for Tutte(x, y) is a ran- 
domised algorithm that takes as input a graph G and a constant e G (0, 1) and outputs a 
value Y such that, with probability at least 3/4, e~ e < Y/T{G; x,y) < e £ . The running 
time of the algorithm is bounded from above by a polynomial in n (the number of vertices 
of G) and e~ x . An FPRAS for PlanarTutte(x, y) is defined similarly. See [7j for further 
details on fully polynomial randomised approximation schemes. 

In earlier work [J], we considered the problem of determining for which points (x,y) 
there is an FPRAS for Tutte (x, y). Our results are summarized in Figure [TJ In particular, 
under the assumption RP ^ NP, we showed the following. 

(1) If x < —1 and (x,y) is not on H or Hi, then there is no FPRAS at (x,y). 

(2) If y < — 1 and (x,y) is not on Hi or H 2 , then there is no FPRAS at (x,y). 

(3) If (x, y) is on H 2 and y < —1 then approximating T(G; x, y) is equivalent in difficulty 
to approximately counting perfect matchings (resolving the complexity of this is a 
well-known and interesting open problem). 

(4) If (x,y) is not on Hi and is in the vicinity of the origin in the sense that \x\ < 1 and 
\y\ < 1 and is in the triangle given by y < — 1 — 2x then there is no FPRAS. 
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Figure 1: The result from [3]. Green points are FPRASable, red points are equivalent to 
counting perfect matchings and gray points are not FPRASable unless RP = NP. We don't 
know about white points. The black line, which is the portion of the hyperbola q = 4 lying 
in y G (—1, 0), is #P-hard. The black points are harder than gray in a complexity-theoretic 
sense. The black region is presumably more extensive than shown. 

(5) If (x, y) is not on Hi and is in the vicinity of the origin and is in the triangle given by 
x < —1 — 2y then there is no FPRAS. 

(6) The two previous intractability results (results (4) and (5)) can be partially extended 
to the boundary of the triangles. 

(7) If (x,y) is in the vicinity of the origin and q = (x — l)(y — 1) > 1.5 then there is no 
FPRAS (excluding the special points at which exact computation is possible). 

An interesting consequence of these results is that, under the assumption RP ^ NP, 
there is no FPRAS at the point (x, y) = (0,1 — A) when A > 2 is a positive integer. Thus, 
there is no FPRAS for counting nowhere-zero A flows for A > 2. This is interesting since 
the corresponding decision problem is in P for example for A = 6. See jl] for details. 

1.4 Approximating the Tutte polynomial of a planar graph 

In this paper we consider the problem of determining for which points (x, y) there is an 
FPRAS for PlanarTutte(x, y). The results of jl] do not help us here because all of the 
constructions are badly non-planar. Our results are summarised in Figure [2} 
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Figure 2: The complexity of the planar case. The shaded gray regions are shown to be 
intractable in Corollary [5j The lower branch of the q = 3 hyperbola (also depicted in gray) 
is shown to be intractable in lemma [6] As Vertigan has shown [12], it is easy to compute 
the polynomial exactly on the hyperbolas q = 1 and q = 2 and at the two special points 
(1, 1) and (—1, —1). (These are shown in green.) 



In particular, under the assumption RP ^ NP, Corollary [5] and Lemma [6] show that 
there is no FPRAS for PlanarTutte(x, y) in the following cases: 

1. x < 0, y < and q > 5; 

2. x < 1, y < 1 and q = 3; 

3. x > 1, y < —1; 

4. y > 1, x < -1. 

For integer q > 4, the point x = 1 — q, y = is of particular interest. As noted earlier 
T(G; x, y) gives the number of proper g-colourings of G. By the 4-colour theorem, there is 
at least one g-colouring, so the corresponding decision problem is trivial, but it is not clear 
whether there is an FPRAS. For q > 5, our result shows that there is no FPRAS for any 
nearby point x = 1 — q/(l+e), y = — e on the hyperbola H q (for any e > 0). However, the 
case of colourings itself (corresonding to the limit point e — 0) remains open. The same 
intriguing situation occurs with the flow polynomial points x = 0, y = 1 — q. 

In a recent posting on ArXiv, Kuperberg (8] independently offers a proof sketch, based 
on the complexity theory of quantum computation, of a result closely related to ours. If 
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the details in the proof sketch can be filled in, then it will strengthen our result in the 
negative quadrant by (i) relaxing the condition q > 5 to q > 4, and (ii) strengthening the 
conclusion to #P-hardness. 

1.5 The multivariate formulation of the Tutte polynomial 

As in jl], we need the multivariate formulation of the Tutte polynomial in order to prove our 
results. The multivariate formulation is also known as the random cluster model [T3l HUj. 
For q G Q and a graph G = (V, E) with edge weights w : E — > Q, the multivariate Tutte 
polynomial of G is defined by Z(G; q, w) = J2AcE w (A)q K ( v ' A \ where w(A) = Y[ e &A w ( e )- 
Suppose (x,y) G Q 2 and q = [x — l)(y — 1). For a graph G = (V, E), let w : E — > Q be 
the constant function which maps every edge to the value y — 1. Then (see, for example 
m (2-26)]) 

T(G; x, y) = (y- iy n {x - 1)~^Z(G; q, w). (2) 

So approximating T(G;x,y) is equivalent in difficulty to approximating Z(G;q,w) for 
the constant function w(e) — y — 1. However, the multivariate formulation is more general, 
because we can assign different weights to different edges of G. 

Consider the following computational problem, which is a planar version of one that 
we considered in [3j. 

Name. MuLTlTuTTE(g; a%, a 2 , a 3 ). 

Instance. A planar graph G = (V, E) with edge labelling w : E — > {«i, «2, ct^}. 
Output. Z(G;q,w). 

Our main tool in proving inapproximability (Lemma [3] below) is showing that 

MuLTlTuTTE(g; a x , a 2 , a 3 ) 

is difficult to approximate if one of ai,a2,«3 is less than —1, one of them lies strictly 
between —2 and 0, and one of them lies strictly outside of this range. 

1.6 Implementing new edge weights 

Let W be a set of edge weights (for example, W might contain the edge weights a±, a 2 
and a 3 from above) and fix a value q. Let w* be a weight (which may not be in W) which 
we want to "implement". Suppose that there is a planar graph T, with distinguished 
vertices s and t on the outer face, and a weight function w : E(T) — > W such that 

w* = qZ st {T)/Z slt {T), (3) 

where Z st (T) denotes the contribution to Z(T;q, w) arising from edge-sets A in which s 
and t are in the same component. That is, Z st (T) = ^ A w{A)q K ^ v,A \ where the sum is 
over subsets A C E(T) in which s and t are in the same component. Similarly, Z s \ t denotes 
the contribution to Z(T;q,w) arising from edge-sets A in which s and t are in different 
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components. In this case, we say that T and w implement w* (or even that W implements 
w*). 

The purpose of "implementing" edge weights is this. Let G be a graph with edge- 
weight function w. Let / be some edge of G with edge weight w(f) = w*. Suppose that W 
implements w*. Let T be a planar graph with distinguished vertices s and t with a weight 
function w satisfying ([3]). Construct the weighted graph G' by replacing edge / with a 
copy of T (identify s with either endpoint of / (it doesn't matter which one) and identify t 
with the other endpoint of / and remove edge /). Let the weight function w' of G' inherit 
weights from w and w (so w'(e) = w(e) if e G E(T) and w'(e) = w(e) otherwise). Then 
the definition of the multivariate Tutte polynomial gives 

Z(G'-q,w') = ^p-Z(G-q,w). (4) 

So, as long as q ^ and Z s \ t (T) is easy to evaluate, evaluating the multivariate Tutte poly- 
nomial of G' with weight function w' is essentially the same as evaluating the multivariate 
Tutte polynomial of G with weight function w. 

Two especially useful implementations are series and parallel compositions. These are 
explained in detail in (5J Section 2.3]. So we will be brief here. Parallel composition is the 
case in which T consists of two parallel edges t\ and e 2 with endpoints s and t and w(ei) = 
Wi and w(e 2 ) = w 2 . It is easily checked from Equation ^ that w* = (1 + Wi)(l + w 2 ) — 1. 
Also, the extra factor in Equation rt2j) cancels, so in this case Z[G'\ q,w') = Z(G; q,w). 

Series composition is the case in which T is a length-2 path from s to t consisting of 
edges e\ and e 2 with w(ei) = wi and w{e 2 ) = w 2 . It is easily checked from Equation ^ 
that w* = wiw 2 /(q + W\ + w 2 ). Also, the extra factor in Equation Q is q + W\ + w 2 , so in 
this case Z(G'; q, w') = (q + W\ + w 2 )Z(G; q, w). It is helpful to note that w* satisfies 

i + X) = fi + J_ > ) ( 1 + JL 

u r J \ W\)\ w 2 

We say that there is a "shift" from (q, a) to (q, a') if there is an implementation of a' 
consisting of some T and w : E(T) — > W where W is the singleton set W = {a}. This 
is the same notion of "shift" that we used in [3]. Taking y = a + 1 and y' = a' + 1 and 
defining x and x' by q = (x — l)(y — 1) = (x' — l)(y' — 1) we equivalently refer to this as a 
shift from (x,y) to (x',y'). 

Thus, the /c-thickening of Jaeger, Vertigan and Welsh [6] is the parallel composition of 
k edges of weight a. It implements a' = (1 + a) k — 1 and is a shift from (x, y) to (x', y') 
where y' = y k . Similarly, the fc-stretch is the series composition of k edges of weight a. It 
implements an a 1 satisfying 

9 (i , 9' ' 



1+ V 

oi \ 01 



It is a shift from (x,y) to (x',y') where x' = x k . 

Since it is useful to switch freely between (q, a) coordinates and (x, y) coordinates 
we also refer to the implementation in Equation ^ as an implementation of the point 



7 



(x,y) = (q/w* + l,w* + 1) using the points 

{(x,y) = (q/w + l,w + 1) | w e 

2 Proving inapproximability 

Our starting point is the following problem. 

Name. Planar cubic Maximum Independent Set. 
Instance. A cubic planar graph G and a positive integer K. 
Question. Does G contain an independent set of size at least Kl 

Lemma 1. Planar cubic Maximum Independent Set is NP -complete. 

Proof. This problem is essentially the same as "Node cover in planar graphs with maxi- 
mum degree 3", which was shown to be NP-complete by Garey and Johnson [21 Lemma 1]. 
Firstly, the complement of a minimum node (or vertex) cover in a graph is a maximum 
independent set. Thus Garey and Johnson's problem is the same as "Maximum indepen- 
dent set in a planar graph with maximum degree 3" . So we just need to show that we can 
transform a planar graph with maximum degree 3 into a cubic graph in such a way that 
the size of a maximum independent set changes in a controlled way. 

It is easily checked that there is a (unique) planar graph T with degree sequence 
(1,3,3,3,3,3). Denote by r the unique vertex of degree 1. Given a planar graph with 
maximum degree 3, we can form a planar cubic graph by attaching (via vertex r) the 
appropriate number of copies of T to the deficient vertices. It is easily checked that each 
copy of T increases the size of a maximum independent set by 2. □ 

We will use the following variant of Planar cubic Maximum Independent Set. 
This variant will help us to maintain planarity in our constructions. 

Name. Planar stretched cubic Maximum Independent Set. 

Instance. A graph G which is the 3-stretch of a cubic planar graph H and a positive 
integer K. 

Question. Does G contain an independent set of size at least Kl 

Lemma 2. Planar stretched cubic Maximum Independent Set is NP-complete. 

Proof. Let m' be the number of edges of H. We claim that the size of a maximum inde- 
pendent set of G is equal to m' plus the size of a maximum independent set of H . 

First, suppose that H has an independent set of size k. We construct an independent 
set of size m' + k in G: For every IN-OUT edge of H the corresponding configuration in G 
can be IN-OUT-IN-OUT. For every OUT-OUT edge of H the corresponding configuration 
of G can be OUT-IN-OUT-OUT. 
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Next, suppose that G has an independent set of size m! + k! for some k! > 0. We 
construct an independent set of size k' in H. Let R be the set of vertices from H in 
the independent set of G and let r = \R\. The independent set of G induces one of 
four possibilities on each edge of H: OUT-OUT, OUT-IN, IN-OUT or IN-IN. Since the 
independent set of G is maximum, all of these four possibilities apart from IN-IN correspond 
to situations in which there is an additional IN vertex in the stretched edge of H. Thus, 
since the size of the independent set is m! + k' we find that m! + k! = r + m! — j, where j 
is the number of IN-IN edges that R induces in H. 

Now take the set R and remove an endpoint from each of these j edges, and we end up 
with an independent set of H of size r — j — k'. □ 

We can now state, and prove, our main lemma. 

Lemma 3. Suppose q ^ [0,5] and that ol\ [—2,0], a 2 G (—2,0) and a 3 < —1. (a^ might 
be equal to a.\ or a 2 .) Then there is no FPRAS for MULTlTUTTE(g; a±, a 2 , CK3) unless 
RP = NP. 

Proof. Suppose H is a cubic planar graph and G is the 3-stretch of H. Let n = \V(G)\ and 
m = \E(G)\, and note that m = |n. Suppose G and K are inputs to Planar stretched 
cubic Maximum Independent Set. We may assume that K < |n, since this is an easy 
upper bound on the size of any independent set in G. 

Our ultimate goal is to construct a planar instance (G', w') of MuLTlTuTTE(g; cti, a 2 , 0:3) 
such that a close approximation to Z(G'; q,w') enables us to determine whether G has an 
independent set of size K. 

As an intermediate goal, we'll construct a weighted planar graph (G = (V,E),w) 
such that a close approximation to Z(G; q,w) enables us to determine whether G has an 
independent set of size K, where w : E — > {(3, a, b} for some conveniently-chosen values 
/3, a and b. Then we'll show how to implement (3, a and b using weights a\, a 2 and 0:3. 
Applying these implementations to the weight graph G will give us G' and w' so that 
Z(G] q, w) is an easily computable multiple of Z(G f , q, w'). 

Let f(x) be the function f(x) = x 3 + 3x 2 . Now in each of the cases that we consider 
(both q > 5 and q < 0), we will define some constants (which depend on q): x, 77, A~ , A + , 
B~ , and B + . The reason for defining them is that we can use these constants (essentially, 
upper and lower bounds on the absolute values of the edge-weights that we will use in our 
gadgets) in order to give a unified proof for the two cases. 

Case 1: q>5. \ = min(l, (q - 5)/6). r] = 3/4. A~ = 1/2. A+ = q. B~ = q. B+ = I0q 3 . 

Case 2: q < 0. x = min(l, \q\). To define the other constants, it helps to make a few 
observations. Let g(y) = /(— 3 — y) and note that g(0) = and that g'(y) < for y > 
so g(y) decreases as y increases from 0. Now let rj > be the real solution of g(rj) = q/2. 
Let A' =3 + i]. Then let y* > be the real solution of g(y*) = q. Let A+ = 3 + y*. Let 
B~ = \q\/3 and let B+ = 4|g|/3 + 2. 

Now, we give quite a few definitions which we will need. Some of them are just using 
the above quantities to form new quantities which we will use as upper and lower bounds 
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on things. Also, we now introduce 5 and s, which may depend upon the input, as well as 
on the above quantities, and are very small. 

Let A* = 1 + 3(A+) 4 + 9(A+) 3 + 3(A + ) 2 + 3A+(1 + \q\) and Q = max(l, \q\). Let 
H = q 2 A+(B + ) 2 , r = \q\{A+) 2 {A + + 3)(B+) 3 , M = max(l, /i,r), and v = 3n - to - 2K. 

Then, using the above definitions, and n and to, we define e as follows. 



e = 



B )' ^,^2-(n+2m+4) Q-Zn+m 



Let L = \q 3 \r]e/2 and R = (B )' /3e. Note that, for sufficiently large n, we have L < 1 
and i? > 1. Define 

^ _ -k y x |y| 

~~ 16 A*3"2 6n M n 2 2m Q 6 ™ ' 

Note that previous estimates on to and K entail v > |n > 1. Now, in both cases, we have 
< rj and 0<5<£<x<l. The definition of 5 ensures that 5 < erj/QA*, which we will 
use later. We will show the following. 

• We can implement a value (5 with |1 + (5\ < 5 (using a±, a<i and 0:3). 

• We can implement a value a with A~ < \a\ < A + and q + e < f(a) < q + 2e and 
|a 2 (a + 3)| > rj. 

• We can implement a value b with B~ < \b\ < B + and — c — 5 < b < — c + 5, where 
c = a 2 + 3a + q. 

The constraints on a imply 1 < A*, \3a + 3a 2 \ < A* and |9a 3 + 3a 4 - 3aq\ < A*. 

Now let Y be a weighted graph with weight function w defined as follows. Y will have 
vertex set V(Y) = {0, 1,2,0,1,2}. The edge set E(Y) of Y consists of three edges (0,0), 
(1, 1) and (2, 2) of weight b and three edges (0, 1), (1, 2) and (2, 0) of weight a. For a fixed 
q, let Z Q \i\2 denote the contribution to Z(Y;q,w) arising from edge sets A in which the 
vertices 0, 1 and 2 are in distinct components. Thus, Zo|i|2 = YIa FLe/i w{e)q K ^ v ^ YS> ' A \ 
where the sum is over all subsets A C E(Y) such that 0, 1 and 2 are all in distinct 
components. Similarly, let Z \ 12 denote the contribution to Z(Y;q, w) arising from edge 
sets A in which the vertex is in one component and the vertices 1 and 2 are in another, 
distinct, component. Finally, let Z i 2 denote the contribution to Z(Y;q, w) denote the 
contribution to Z(Y; q, w) arising from edge sets A in which the vertices 0, 1 and 2 are all 
in the same component. 

Let 



n 2 7 K-l 

Q ^012 



Zq\1\2 



-K l z o|i|2| \q\ X ■ 



From the definition of Zo|i2 we see that 

Z \i2 = q 2 ab 2 (c + b) 
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so \Z \ 12 \ < 5ji. Similarly, 

Zo\x\2 = Q 3 ^ + 36 2 (2a + q) + (36 + g)(a 3 + 3a 2 + 3ag + q 2 )) . 
so for some d G [—5, S] and e G [e, 2e], explicitly, e = a 3 + 3a 2 — g and d = a 2 + 3a + q + b, 



Z 0W 2 = -gV(a + 3)(a 3 + 3a 2 - q) + d 3 g 3 + d 2 (-3a - 3a 2 )g 3 + da 3 (9a 3 + 3a 4 - 3ag). 

= -g 3 a 2 (a + 3)e + d 3 g 3 - d 2 (3a + 3a 2 )g 3 + dg 3 (9a 3 + 3a 4 - 3ag). 

□ o, „, / d 3 d 2 (3a + 3a 2 ) d(9a 3 + 3a 4 - 3ag) 

= -g 3 a 2 a + 3 e - + > - 2 , 

\ a^(a + 3) a^(a + 3) a z (a + 3) 



Now, the absolute value of each of the right-most three terms is at most 5 A* ft] and by 
the definition of 5, this is at most e/6. Thus, 

L < |g 3 a 2 (a + 3)| E - < \Z 0]ll2 \ < \q 3 a 2 (a + 3)|3e. (5) 

Also Z i2 = ga 2 (a + 3)6 3 so \Z i 2 \ < r and 

|Z 012 | > |ga 2 (a + 3)|fi- 3 . 



Finally, 

q 2 \Zoi2\ q 2 \qa 2 (a + 3)\B' 3 = 

\Z0\1\2\ ~ |g 3 a 2 (a + 3)|3e ' 1 ] 

Here is the construction of G = (V,E). See Figure [3J First, fix any ordering on the 
vertices of H. 

Next, let's set up some useful notation for the graph Y, which we will use as a gad- 
get. A particular copy Y x of this gadget will have vertex set V x = {(x,0), (x, 1), {x,2}, 
(x, 0), (x, 1), (x, 2)} with vertices (x, 0), (x, 1) and (x, 2) arranged in clockwise order around 
the outer face. The edge set E x consists of three edges ((x, 0), (x, 0)) , ((x, 1), (x, 1)), and 
((2, 2), (x, 2)) of weight b and three edges ((x, 0), (x, 1)), ((x, 1), (x, 2)) and ((x, 2), (x, 0)) 
of weight a. We will construct a copy Y u of the F-gadget for every vertex u G V(H). 

Now, associate every edge (u,v) of H with two indices i UtV and in {0,1,2} in such 
a way that the graph with vertex set {J ueV (H) ^ U anc ^ e dge set 

|J E U U |J «u, <u,z»,„)) 

is planar. We will construct two copies, and Y vu , of the F-gadget for every edge (u, v) 
of H. These correspond to the vertices of G along the two-stretched edge (u, v) of H. 
Thus, we have one Y gadget for every vertex of G. 
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Figure 3: The portion of G corresponding to edge (u, v) of H. In the picture, we assume 
i U)V = and i V)U = 1. Edges of E are depicted as solid black lines and edges of E' are 
depicted as dashed red lines. Where two vertices have been identified, both the original 
labels have been displayed. 

The vertex set V is constructed from Uuev(-ff) V U ^V}( U v)&E(H)(y UV ^V VU ) by identifying 
some vertices, In particular, for every edge (u,v) of E(H) with u < v , identify (u, i u>v ) 
with (uv,0). Also, identify (uv,2) with (vu, 1). Finally, identify (vu,0) with (v,i V)U ). Note 
that G has one vertex for each vertex of H and two vertices for each edge of H so n = 
\V(H)\+2\E(H)\. Also m = 3\E(H)\. So \V\ = 6\V(H)\ + 12\E(H)\-3\E(H)\ = Qn-m. 

Let E = [J U&V{H) E u U [j(u,v)€E(H)( EUV u Em )- E is a11 of the internal edges in the Y 
gadgets. So \E\ = 6\V(H) \ + 12\E(H)\ = 6n. Let E' be the set of m edges with weight (5 
constructed as follows. For each edge (u, v) of E(H) with u < v, let i — i U)V — 1 mod 3 and 
let j = i VtU — 1 mod 3. Add edges ((u, i), (uv, 1)), ((uv,l), (vu,2)) and ((vu,l), (v,j)) to 
£". Let E = EU E'. Note that G is planar. Also, 

Z(G;q,w) = J2J2 w(A)w(B)q K ^' AuB \ 

ACE BCE> 

Let (S,D,T) be a partition of V(G) into three sets 5, D and T, some of which could 
be empty. (In the following, it will help to think of S as "singleton", D as "doubleton" 
and T as "triple".) Let As,d,t denote the set of subsets ACE such that the following 
statements are true in the graph (V, A). 

• For every x G S, the vertices (x, 0), (x, 1) and (x,2) are in a single component, 

• for every x £ D, the vertices (x, 0), (sc, 1) and (x, 2) are in two distinct components, 
and 
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• for every x G T, the vertices (x, 0), (x, 1) and (x, 2) are in three distinct components. 
Let Zs,d,t be the contribution to Z(G; q,w) from A G As^d.t- That is, 

^w= E E ^w^^), 

so Z(G; q, w) = ^2(g D T) Z.s,d,t> where (S, D, T) ranges over all ordered partitions of V(G). 
First, fix a partition (S,D,T) in which D is non-empty. Then 



f)6n 

\z s , d ,t\< E l^&EM 5 )! 

and 

E |™(A)<f ^| < 2 6 "|Z | 12 | IZqizI"" 1 < 2 6n ^r"- 1 . 

Now, each of the 2 m choices for B contributes at most \w(B)\ < \(3\ m < (1 + <5) m < 2 m . 
So, using the definition of S and the lower bound in Equation ([5]), 

E l^wl < T2 6n 5iir n - l ^-2 2m < M//16. (7) 

Next, fix a partition (S,D,T) in which .D is empty and S is not an independent set 
of G. We give an upper bound for \Zs,d,t\- Suppose (x,y) is an edge of G with x G S 
and y £ S. Let e be the edge between {(x, 0), (sc, 1), (x, 2)} and {(?/, 0), (?/, 1), (y, 2)} 
that is in E' . For any B C £", let B' = B — {e}. Note that, for any set A G As,d,t, 
k(V",4 U 5') = U B' U {e}). To see this, note that (x,0), (x, 1) and (x,2) are 

in the same component of (V, A) (since x G S) and the same is true of (y, 0), (y, 1) 
and (y, 2); moreover, some vertex in {(sc, 0), (x, 1), (x, 2)} is identified with some vertex in 
{(y,0),(y,l),(y,2)}. Thus, 

Zs,d,t — E E w (A)w(B)q K ^> AuB) 

AeA s ,D,T BCE' 

= {i + P) E E M^MW^ , 

and hence \Zs,d,t\ is at most 

+ E l^^lrk E K5 / )l<^2 6 v^2 2 - 

So, using the definition of S and the lower bound in Equation ([5]), 

E I Z s,At| < 2 n 52 6n T n f-^2 2m < tf/16. (8) 



D=0 
5 not independent 
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Finally, fix a partition (S, D, T) in which D is empty and S is an independent set of G 
of size k. 

Z w = w{A)q< 9 ^ w(B)q^^uB)- K( v,A) 

AeA S ,d,t BCE' 

„-myk yn-k \^ n\B\ k(V ,AUB)-k(V ,A) 

— V ^012^01112 2-^1 p ** 
BCE' 

Now let V = V - \J xeV{G) {{x,0), (x,l), (x,2)} so \V'\ = \V\ - 3n = 3n - m. Let 
T = (V',E') and let T s be the graph constructed from T by identifying (x,0), (x, 1) and 
(x,2) if x G S. Note that |V(r s )| = 3n-m-2k. Then for any B C £7' and any A G A s ,d,t, 
k(V, A) - k(V, AUB) = k(V(T s ), 0) - k(V(T s ), B). We conclude that 



Zs,d,t - 3 m ZQ l2 Z^ 2 /5 |jB| g K 



i7 «(v(r s ),i?)-«(v(r s ),0) 

BCE' 

_ o-m 7 A: yn-fc -3n+m+2fc n\B\ k(V{V s ),B) 

— O Zj 012 Zj Q ^ 2 q ^ P H 

BCE' 

Then we have ^ 

Z S ,n,T = (^) Z^ 2 q~ 3n Z(T s ;q,w), (9) 

where w is taken to be the restriction to E', i.e., the constant (3. 
Now we want some bounds on 

Z(T s] q,w) = J2p\B\ q «V(r s ),B)_ 

BCE' 

A crude upper bound is \Z(T s ;q,w)\ < 2 m (l + 5) m Q 3n - m < 2 2m Q 3n ~ m . So, using the 
definition of e, 



^ \Zs,d,t\ < 2 r 



D=0, |S|<if, 
S is independent 



q 2 Z 



012 



^0|1|2 



K-l 

n LI" 3Ti o2mr)3n 



^oii| 2 rigr 2 2m Q jn - m < */i6 (io) 



For the lower bounds, we fix k = K so |V(rs)| = v. 

Case 1: q > 5. For a lower bound on \Z(Ts;q,w)\ and information about its sign, we 
use [T4"t Theorem 1]. (The method can be traced back to pQ.) Since x — (<? ~~ 5)/6 and 
5 < x, we have g > 5(1 + 5) + x > 5|/3| + Now the inductive proof of [TJJ Theorem 1] 
gives (for some vertex t> of of degree at most 5) 

Z(T S ; q, w)>(q- 5\P\)Z(T S - v; q, w) > (q - 5\(3\f iFs)l > X u - 

Thus, for q > 5 we have shown that the sign of Z(T,s',q, w) is the same for all S (it is 
always positive). Also, we have shown 

\Z(T s ;q,w)\>x' / . (11) 
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Case 2: q < 0. To determine the same facts for q < we use Theorem 4.1]. Note 
that T s has no loops. Let Ci,...,C u denote the coefficients of Z(T s ;q,w), viewed as a 
polynomial in q, so Z(T S ; q, w) = Y^j=\ ^0 ' ■ Let ttj = 1 if Cj > 0, ttj — if Cj = and 



7Tj = — 1 if Cj < 0. Then 



z(r s ;q,w) = (-irj2(- i r j ^\ c . 

3=1 



Jackson and Sokal [3, Theorem 4.1] showed (assuming 5 < 1, which implies —2 < /? < 0) 
that {-l) u -Hj > 0. So 



z(r 5 ; q,w) = (-if £ 1^1 
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i j 



Note that for j = z/, Cj = 1 so Equation (11 ) holds and the sign of Z(T S ; q, w) is the same 



for all size-i^ sets S (the sign depends on the parity of v). 

Now, for a partition (S,D,T) in which D is empty and S is an independent set of G 
of size K, Equations (|9]), ^ and (11) give 



J S,D,T 



> (12) 



Since the sign of ^(r^; q, w) is the same for all S, it is apparent from ^ that the sign of 
Zs,d,t depends only on the sign of q, the sign of Z 012 , the sign of Z \i\ 2 and the parity of 
K and n. It does not depend on the set S. 

So if G has iV > independent sets of size K then by Equations Q, @, (JlO]) and ([l2|), 



\Z(G;q, w)\ > — 3\l//16 > 3^/4. On the other hand, if G has no independent sets 
of size K then the same equations give \Z(G;q, w)\ < 3^/16 < \&/4. So if we could 
approximate Z(G; q, w) within a factor of | then we could determine whether or not G has 
an independent set of size K. 

To finish the proof, we must first show how to implement /3, a and b using a\, a 2 
and a>3. The first of these is already done in [4]. To implement /3, choose a positive integer 
k such that \{a 2 + l) k \ < 5 then implement j3 by fc-thickening a 2 - Before showing how to 
implement a and h it helps to have a look at the implementations that we can obtain by 
series and parallel compositions. 



2.1 Moving around hyperbolas: Case 1: q > 5 

The weights that we have available for our implementations are a%, a 2 and 013. It will 
be useful to use (x,y) coordinates as well as (q,a) coordinates since series compositions 
power x and parallel compositions power y. Thus, we define y' ri = 1 + a>i and x\ = q/{y[ — 
1) + 1 for i E {1, 2, 3}. (The primes are just there because we use yi for something else.) 

The main purpose of the section will be to show how to use the values y[, y' 2 and y' 3 
to implement a value that is very close (within ±ir for some exponentially small value tt) 
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to some target value T > 1. To do this efficiently (without exponentially big gadgets) we 
need to be somewhat careful about decomposing T. Note that we won't actually use y' 2 in 
this implementation. 

Let (xi, Di) be the point on the hyperbola (x\ — — 1) = q given by y\ = y[ 2 . Note 
that yi > 1 so X\ > 1 and that we we can implement (x\, yi) by 2-thickening from (x[, y[). 
Let 

q 



q. Note that, for every integer 



X X 3 ~ 1 + ' 

Let Xj be the corresponding value so that (xj — l)(yj - 
j > 1, we can implement 
% > 1 and % > 

Now, for every integer j > 1, we recursively define a quantity <ij in terms of the values 
of g?i, . . . , In particular, 



Xj, yj) by j-stretching from (x\, y\). Also, since x\ > 1, we have 



(I; 



iog(rn£}w- 



log(% 



Let = rifci Note that T/y m < y'^ < T. Also, since c?j is a non-negative integer, 
we can implement y!^ by ^-thickening ye and then combining these in parallel. 

First, here is how to implement a value y which is within ±7r of some value T > 1 for 
some < Ti < 1: Let 

"log(gT/7r + 1) 



m 



log(xi) 
vr/T), so l/y m 



> 1 — 7r/T. Let y = y'm and note that 



Note that y m < 1 + vr/T < 1/(1 
T - 7T < y < T. 

To see that this implementation is feasible, note that m is not too large. In particular, 
for fixed q and y[ and upper and lower bounds on T, m is bounded from above by a poly- 
nomial in the logarithm of 7r _1 . To finish, we must show that the same is true of d\, . . . , d m . 
Here, the key observation is that y^ < T/y'j_ 1 < yj-%, so dj < log(yj_i)/log(?/j). Then for 
Vj < 5/4, say, we have |(?/j — 1) < log(yj) < yj — 1 which suffices. 

Second, here is how to implement a value y which is within ±7r of some value T < — 1 
for some < tt < 1: Use the above method to implement a value y' satisfying 



|T| 



7T 



103 1 



<y'< 



|T| 

551 



Then take y' in parallel with j/3 to get a value |/ satisfying T < y < T + tt. 



2.2 Moving around hyperbolas: Case 2: q < 



This is very similar to Section 2.1, Here, we start by implementing a point (xi,yi) with 
Xi < —1. Then we just use odd values of j and it suffices to take 

"log(|g|T/7r)" 



m 



log(|xi|) 
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The point (xx,yi) is reached as follows. If y' 2 < 1 + |g|/2 then we can take (xi,yi) 
(q/(y' 2 — 1) + l,y'i 2 ) since X\ < —1. Otherwise, proceed as follows. Let 



2 l + |?|/2 



Choose a positive integer j so that 



- 1) + 1) ; 



< 0. 



There is such a j since y 2 £ (~ 1> !)• Now let (£, y) be the j-stretch of (x^y^) so 1 — £ < 
y < 1. Now let 

_ log((l + |g|/2)M 2 )' 



log(y) 



Note that A: is a positive integer since y' 2 > l + |g|/2. Let (xx, y%) be the parallel composition 
of (q/{y[ 2 — 1) + l,y'i 2 ) with the /c-thickening of (x, y). Thus, y\ = y k y' 2 ■ Note that 
1 < y(l + |g|/2) < yi < 1 + |g|/2 so si < -1. 



2.3 Implementing a: Case 1: g > 5 

We start by noting that any a that satisfies 

q + e < f(a) <q + 2e 



(13) 



also has \ < a < q so A < \a\ < A + and |a 2 (a + 3)| > rj. So we just need to see how to 
implement a value of a that satisfies (13). Suppose that for some positive rational number 



a', a' is too small to be a suitable value of a. In that case, we will argue that a 1 + e 2 is not 
too large to be a suitable value (hence, there is a suitable value in an interval of width e 2 ). 
Thus, we can finish by using the argument in Section 2.1 taking it = e 2 and T to be the 
solution to the equation f(T) = q + 2e which we know satisfies A~ < T < A + . 

To see that a' + e 2 is not too large, suppose f(a') < q + e so (since q > 1) we have 
a' < q. Then f{a' + e 2 ) - f(a') = (6a' + 3a' 2 )e 2 + (3 + 3a> 4 + e Q . Now using the definition 
of e (and assuming that n and m are sufficiently large) we see that this is at most e. 



2.4 Implementing b: Case 1: q > 5 

First note that B~ < c — 1 < c — 5 and c + 5<c + l< B + so it suffices to implement 
a value b with — c — 5 < b < — c + 5. Once again, we use the argument in Section |2.1 
choosing T = c and n = 5. 
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2.5 Implementing a: Case 2: q < 

Start by noting that if y satisfies /(— 3 — y) < q + 2e < q/2 then y > rj. Also, if y satisfies 



/(— 3 — y) > q + e > q then y < y*. Thus, if a = — 3 — y satisfies (13) then A < \a\ < A + 
and \a 2 (a + 3)| > i]. So we just need to see how to find a value of a that satisfies (13). 



This is now essentially the same as Section 2.3 



2.6 Implementing b: Case 2: q < 

By equation (p~3 

n 4- If 

<a 2 + 3a< 

a 

so since a < —3, 



q + 2e^_ 2in _^q + e 



,g + 1 < +q < 1 +q < a z +3a+q < +q < q 1 + - < -- \q \< 0. 

6 J a a a \ a J 3 

Thus, —c — 5 > B~ and — c + 5 < B + . So it suffices to implement a value 6 with 



— c — 5 < b < — c + 5. For this, we just use the argument in Section |2.2[ choosing T = c 
and 7r = 5. 

□ 

Theorem 4. Suppose (x,y) G Q 2 satisfies q — (x — l)(y — 1) G" [0,5] . Suppose also that 
it is possible to shift the point (x,y) to a point (xi,yi) with y\ G" [—1,1] and to a point 
(22,2/2) with y2 G (—1, 1) and to a point (23,2/3) with y$ < 0. JTien there is no FPRAS for 
Tutte(x,?/) unless RP = NP. 



Proof. This follows easily from Lemma [3] and is similar to the proof of [4, Theorem 2]. For 
completeness, here is a proof. 

Let a = y — 1 and a* = yi — 1. Let (Tj, e*) be a graph that shifts (x, y) to (xj, yj). Note 
that (Tj,ej) shifts (q,a) to (g, aj). 

Suppose (G, w) is an instance of MuLTlTuTTE(g; ai, 0:2, 03) and note that ai, «2 and «3 
satisfy the conditions of Lemma [3} Suppose that G has mj edges with weight c^. Denote 
by G the graph derived from G by applying the above shifts — replacing each edge with 
weight ctj with a copy of Tj. Let w be the constant weight function which assigns weight 
a to every edge in G. Then by Equation [IJ 

Z(G;q,w)=\^ ^ 1 g2 j g2 j ^(G;?,uO, 



so by Equation [2j 



v,/. , / ^s|t(Ti)\ mi /Z s | t (T 2 )\ m VZ s | t (T 3 )\ m3 ''''/■/ / " ' 

Z(G;g,w) = ^ ^ ^ (y- 1) (x - 1) T(G;x,y), 



q 2 J \ q 2 J \ q 2 



where n is the number of vertices in G, and k is the number of connected components in 
G. 
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Thus an FPRAS for Tutte(x, y) would yield an FPRAS for MuLTlTuTTE(g; ai, a 2 , a 3 ), 
contrary to Lemma [3} 

□ 

The following corollary identifies regions where approximating Tutte(x, y) is intractable. 
It is illustrated in Figure [2} 

Corollary 5. Suppose RP ^ NP. Then there is no FPRAS for PlanarTutte(x, y) when 
(x,y) is a point in the following regions, where q denotes (x — l)(y — 1): 

1. x < ; y < and q > 5; 

2. x > 1, y < -1; 

3. y > 1, x < -1. 

Proof. We will show that for each point (x, y) in the following regions, we can shift to 
points (x\,yi), (22,2/2) and (23,1/3) satisfying the conditions of Theorem [I| 

For the remaining cases, we use the fact that, when G is a planar graph and G* is any 
plane dual of G, T(G;x,y) = T(G*;y,x) [T31 §3.3.7] (so the fact that there is no FPRAS 
at (x,y) implies that there is no FPRAS at (y,x) and vice- versa). The regions that we 
consider are as follows. 

1. x < —1, y < — 1 and q > 5: We can take (23,1/3) and (21,3/1) to be (x,y) since 
y < —1. We can realise (22,^2) using a large, odd, /c-stretch so y 2 = q/(x k — 1) + 1 
which is in the range (—1,1). 

2. — 1 < 2 < 0, y < —3/2 and q > 5: Note that the condition y < —3/2 is implied by 
the bounds on 2 and q. As above, we can take (23,^3) and (21,1/1) to be (x,y) since 
y < — 1. Next, let T = —5/4 and let tt = 1/8. Use the method from the very end 



of Section 2.1 to shift to a point (x",y") with T < y" < T + tt. This method uses 
the points (21,1/1) and (23,1/3). Now note that, since q = (x" — l)(y" — 1) > 5, both 
2" < —1 and y" < —1 so we can use the previous case to shift to a suitable (22,2/2)- 

3. x > 1, y < —1: Note that q < 0. In this case we can take both (21, y\) and (23, 2/3) 
to be the point (2, y) since that gives y\ — y% < — 1. We get to (22, 2/2) by a j-stretch, 
for sufficiently large j. This gives 

□ 
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3 The lower branch of q = 3 



The following is NP-hard [3J. 

Name. Planar 3-Colouring. 

Instance. A planar graph G. 

Question. Does G have a proper 3-colouring? 

The following lemma gives hardness for approximating the Tutte polynomial on the 
lower branch of the q = 3 hyperbola. See Figure [2} 

Lemma 6. Suppose RP ^ NP. Then there is no FPRAS for Tutte(x,?/) when (x,y) 
satisfies (x — l)(y — 1) = 3 and x,y < 1. 

Proof. We will consider a point (x,?/) with — 1 < y < 1. The remaining cases follow by 
symmetry between x and y as in the proof of Corollary [5j Let G = (V, E) be an input to 
Planar 3-Colouring with n vertices. For an even positive integer k, let G k be the graph 
formed from G by fc-thickening every edge and let E k be its edge set. It is well-known (see, 
for example, jH Section 5.1]) that 

T{G k - x, y) = (y- - l)-*^) £ y~^\ 

<t:V->{1,2,3} 

where mono(cr) is the number of edges in E k that are monochromatic under the map o. 
Note that mono(cr) is an even number, since k is. Thus, J2a-v^{i 2 3} ?/ mono(CT ) is a positive 
number which is at least 1 if G has a proper 3 colouring and is at most 3 n y k otherwise. 
Choosing 

k= f log(4-3") l 
log (l/y) 

we have 3 n y k < 1/4, so a 2- approximation to T{G k ;x,y) would enable us to determine 
whether or not G is 3-colourable. □ 
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